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ABSTRACT 

The gravitational potential of clusters of galaxies acts as a cosmic telescope allowing us to find and 
study galaxies at fainter limits than otherwise possible and thus probe closer to the epoch of formation 
of the first galaxies. We use the Bullet Cluster 1E0657— 56 {z = 0.296) as a case study, because its 
high mass and merging configuration makes it one of the most efficient cosmic telescopes we know. 
We develop a new algorithm to reconstruct the gravitational potential of the Bullet Cluster based on 
a non-uniform adaptive grid, combining strong and weak gravitational lensing data derived from deep 
HST/ ACS F606W-F775W-F850LP and ground-based imaging. We exploit this improved mass map 
to study z ~ 5 — 6 Lyman Break Galaxies (LBGs), which we detect as dropouts. One of the LBGs is 
multiply imaged, providing a geometric confirmation of its high redshift, and is used to further improve 
our mass model. We quantify the uncertainties in the magnification map reconstruction in the intrinsic 
source luminosity, and in the volume surveyed, and show that they are negligible compared to sample 
variance when determining the luminosity function of high-redshift galaxies. With shallower and 
comparable magnitude limits to HUDF and GOODS, the Bullet cluster observations, after correcting 
for magnification, probe deeper into the luminosity function of the high redshift galaxies than GOODS 
and only slightly shallower than HUDF. We conclude that accurately focused cosmic telescopes are 
the most efficient way to sample the bright end of the luminosity function of high redshift galaxies 
and - in case they are multiply imaged - confirm their redshifts. 

Subject headings: cosmology: dark matter - gravitational lensing - galaxies:clusters:individual:lE0657- 
56 - galaxies: high-redshift 
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1. introduction 

Studying galaxies at high redshifts is crucial for un- 
derstanding both their formation and evolution, and the 
role they played in cosmic reioniza tion. However, recent 
observations of z > 6 objects (e.g., iBouwens et al.|[2QQ8f ) 
show that UV-bright galaxies alone are not believed to 
be sufficient to reionize the Universe (unless for exam- 
ple, the escape fractions or clumping factors are dramat- 
ically different from current assumptions). It is therefore 
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important to continue to find radiation sources at the 
highest redshifts and study their physical properties. 

Observations of galaxies at these high redshifts are 
challenging, not only due to the large luminosity dis- 
tance to these objects, but also due to their lower lu- 
minosity, since the luminosity function evolves com- 
pared to galaxies at moderate redshifts ~ 2; see e.g. 
Reddy & Steidel 2009, Bouwens et al. 2007). Recently 
our knowledge of the early cosmic epoch of galaxy forma- 
tion {z > 4) has vastly increased, providing an opportu- 
nity foru^_t^^ (Stark e tld] 
'2009', 'Bouwens et al." '20081 |2QQ7i [Y oshida et al.. .20061 . 
Stanway et al. 2005, Giaval isco et a l. 2004b, and refer- 
ences therein). Comparing the measured luminosity 
function with results from simulations provides insight 
on the efficiency of star formation as a function of mass 
scale, an d on the feedback mechani sms that may limit it 
(see e.g. Ivan den Bosch et al.ll2003D . 

We can find high-redshift galaxies by searching for 
the redshifted Lyman break using broad band photom- 
etry (Steid el et al.lll996f ). These so urces, known as Ly- 
man Break Galaxies (LB Gs, see e.g. IVanzella et al.|[2'009l 
for recent studies, and iGiavaliscol l2002l for a review) 
are the best studied and the largest sample of galax- 
ies at redshifts z ^ 5. One can identify z 5 ob- 
jects by t heir non- detection in th e V-band and blue- 
ward (e.g. iGiavalisco et al]|2004al [Stanw a v et al] I2005L 
iHildebrandt et all 120091 . IStark et al]l2009[ ):"such objects 
are referred to as "F-band dropouts." Similarly, ob- 
jects di t z 6 are associated with i-band non-det ection 
(e.g. IStanwav etHI 12001 IBouwens etHI I2008[ ). and 
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z 7-8 with z-band ( e.^. iBouwens fc IllingwQrthll2QQd 
Bouwens et all I2QQ8L iHenrv et all I2QQ9I . lOesch et all 
2QQ9bl ). To date, preliminary results usi ng the new 
WFC3/ IR d ata have now been released. lOesch et al.l 
([2009a) and Bunker et al. (2009) also study t he z-band 
dropouts, whereas Bouwe ns e t al. (2009t^, lYan et aP 
|2009), and lMcLure et all ([20091 ) study F-band dropouts 
which are potential z ^ 8 — 8.5 objects. At even 
higher redshifts , z ^, S objects are observed as J- band 
dropouts (e.g. lYan et all 120091 . iHenrv et al.1 120081 . and 
iBouwens et al.l [2005f ).^^ The main limitations of these 
experiments to date are the relatively low expected num- 
ber counts (even in the deepest fields), and contamina- 
tion {z ^ 2 dusty galaxies and L and T dwarf stars can 
also share the characteristic dropout appearance). Fur- 
thermore, the depths probed by these surveys are around 
characteristic luminosity for high-redshift population 
and do not reach far into the low luminosity regime where 
the faint-end slope is constrained, especially at higher 
redshifts. 

We argue in this paper that at the luminosity regime 
down to L L^^ determining the high redshift lumi- 
nosity function can be performed more efficiently when 
using galaxy clusters as cosmic telescopes with respect 
to blank fields. This technique was proposed shortly af- 
ter the first grav itationally-len sed arcs in galaxy clusters 
were discovered ([Soucaill ll990) and has been successfully 
applied to ever- improv i ng data (e.g. , Zheng et al. 2009, 
[Bouwens etaB l2009al iRichard et all I200a lEllis et al.1 
|2001'). Observations of galaxy clusters used as cos- 
mic telescope is consistently delivering reco rd holders in 
the s e arch for the hi ghest redshift galaxies ([Kneib et al.l 
I2OOI iBradlev et al.1 12008). The reason for this is that 
the magnifying power of a massive cluster lens typically 
allows the detection of objects more than a magnitude 
fainter than the observation limit. Due to this same mag- 
nification, the effective solid angle of the survey volume 
decreases; however, since the luminosity function is prac- 
tically exponential at the magnitudes we probe, the lens- 
ing magnification gives us a substantial net gain. In- 
deed, if the faint end slope o f the luminosity fun ction 
at high redshift is steep (e.g. IBouwens et al.|[2007f ). we 
gain at magnitudes below and around the characteristic 
magnitude M*. One concern is that we need to prop- 
erly account for the errors introduced when we determine 
the true o bserved volume with a cosmic telescope (e.g. 
IBouwens^ et al. 2009a). However, we quantify the uncer- 
tainties in the magnification map and show that they are 
negligible compared to sample variance when determin- 
ing the luminosity function of high-redshift galaxies. 

Cosmic telescopes offer two further advantages. First, 
the detected sources behind a cluster lens are magnified 
in apparent size. As a result, we can resolve smaller phys- 
ical scales than would otherwise be possible, and begin 
to actually measure the properties of z > 6 galaxies on 
a case-by-case basis. Resolving the galaxies also helps in 
rejecting contamination by foreground cold stars. Sec- 
ond, some of the sources will be multiply imaged. If this 
is the case, then they can be readily distinguished from 

The z ~ 9 candidate object from ' Henry et al'1l2QQ8l was later 
confirmed to be a lower redshift one and the number of candi- 
dates from Bouwens et al. 2005 appears consistent with the ex- 
pected contamination from low-redshift interlopers. 



the main contaminants, typically intermediate redshift 
galaxies and cold stars. The positions, where multiply 
imaged systems form are redshift dependent, so a z > 5 
multiply imaged source can be readily distinguished from 
(multiply-imaged) z ^ 2 source and singly imaged stars. 
Thus, in the magnitude range we observe here, a cos- 
mic telescope survey is more efficient than a competing 
imaging survey in unlensed fields to the same depth. 

The most effective gravitational lenses for such stud- 
ies are massive and highly elongated (in projection, on 
the plane of the sky) galaxy clusters. It is easy to un- 
derstand the first requirement, as more massive objects 
are more powerful lenses. The ellipticity is important as 
well, as the area of highest magnification is located at 
the semi major-axis intersection with the critical curve 
(i.e. curve connecting points of formally infinite mag- 
nification). Equivalently, sources close to the major axis 
cusp have the largest magnific ation. The clu ster of galax- 
ies 1E0657— 56, discovered bv I Tucker et al.l (|1995), is one 
of the hottest, most X-ray luminous c lusters known. It 
is also a plane-of-the sky merger (e.g. iMarkevitch et all 
I2OO2L LBarrena et alj 120021 ) , hence giving a highly elon- 
gated mass distribution in projection. 

It is also at a favorable redshift, which contributes to 
make it a very efficient lens. The critical curve radius 
decreases with increasing lens redshift for a source at a 
fixed (higher) redshift. One might thus think that lower 
redshift clusters 0.05 < z <0.2 are better suited for these 
surveys. However, a small solid angle coverage by cluster 
member galaxies is desirable, hence reducing the effec- 
tiveness of clusters that are at too low redshift. Whereas 
lensing efficiency scales with angular diameter distances 
between observer, source, and the lens, the size of cluster 
members only depends upon the latter. There is an op- 
timal balance between the size of the critical curve and 
the size of the cluster members. In addition, a good cos- 
mic telescope candidate should also show many strongly 
lensed images (i.e. efficient also for lensing of 2: ~ 1 — 3 
sources), which are needed to reliably reconstruct its 
magnification map. Bullet Cluster, with > 10 definite 
strongly lensed systems and the lens redshift of z ^ 0.3, 
is well suited for these studies. Specifically, since we are 
trying to detect objects close to the detection threshold, 
we do not want to rely on perfect image subtraction of 
the cluster members. For clusters at much higher red- 
shifts, on the other hand, the critical curves (and the 
area of high magnification) are not large enough. 

A key requirement for a cosmic telescope survey is an 
absolutely calibrated, high-resolution mass reconstruc- 
tion for the lens. This is necessary to accurately con- 
vert the observed solid angle into the actual cosmic vol- 
ume, as well as to reconstruct the source fiuxes and 
so determine their luminosity. In this paper we will 
describe a newly developed technique to reconstruct 
the mass distributions of galaxy clusters. As we will 
show, the technique allows us to improve the accuracy 
in the very center of this cluster, which is crucial to 
accurately predict the magnifications of dropouts (and 
other objects in the field). The key feature of this 
new technique is the optimal use of information from 
both multiple image systems (strong lensing from deep 
high-resolution HST/ ACS data) and from the distortions 
of the singly-imaged background sources (weak lensing, 
from combined ACS and ground based data) using adap- 
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tive grid. Multi-resolution methods have been used in ei- 
t her strong or weak gravitational lensing in the past (see 
e.g. lM arshaH"2QQ6' 'Diego et al."2QQ7L lMerten et al.ll2QQ8l . 
iDeb et al. 2008, Julio & Kneib 200l). The new method 
described here, uses information from both strong and 
weak lensing regimes, on a multi-level grid that allows 
us to efficiently reconstruct features at the maximum 
resolution allowed by the data. This new method and 
map supersedes our previously-inferred potential map, 
which was limited in its resolution due to the inability of 
the mass-modeling routine to simultaneously reconstruct 
large angular scale features (constrained by the relatively 
low signal-to-noise weak lensing data), and the smaller- 
scale features in the high signal-to-noise strong lensing 
regime. 

This paper is structured as follows. Section [2] describes 
the improved method to reconstruct the mass distribu- 
tion on a non- uniform grid of pixels. Section [3] describes 
the data and its reduction procedures, and in Section [H 
we present the mass reconstruction of Bullet Cluster. 
Section [5] then gives the results of our search for z ~ 5 — 6 
objects using our newly-calibrated lens, and discusses 
further the advantages and disadvantages of performing 
such studies behind clusters. Our conclusions are sum- 
marized in Section [6l 

Throughout the paper we assume a ACDM cosmology 
with = 0.3, flA = 0.7, and Hubble constant Hq = 
70 kms-^ Mpc~^ At the cluster redshift z^i = 0.296, 
the physical length scale is 4.4 kpc/arcsec. All the co- 
ordinates in this paper are given for the epoch J2000.0. 
All magnitudes are given in the AB system. 

2. STRONG AND WEAK LENSING MASS 
RECONSTRUGTION ON A NON-UNIFORM ADAPTED 
GRID 

Our combined strong and weak lensin g analysis fol- 
l ows t he algorithm first proposed by Ba rtelma n n et al ] 
(|1996f ). The impl ementation that inc ludes strong lens- 
ing is described inlBradac et al .' ('2005') and implemented 
on ACS data in lBradac et all (2006. 2008a, b). However, 
in these works the method showed limitations due to the 
fact that the resolution of the reconstruction was uniform 
across the observed field. Therefore, all our reconstruc- 
tions were limited to small fields (a few arcmin on a side). 
By using a reconstruction grid whose pixel scale varies 
across the field, we are able to overcome these limitations 
and map the gravitational potential of a cluster lens over 
a much larger field of view (Bradac et al. 2009, in prepa- 
ration). More importantly for this paper, we are able to 
achieve increased resolution in the cluster centre (close 
to where we see strongly lensed images), and hence the 
magnification map in the regions of high magnification 
is more accurate. 

All lensing observable quantities can be written 
as spatial derivatives of the lens potential (see e.g. 
iBartelmann fc Schneider! l2QQlf ). In the past we have 
modeled a cluster's projected gravitational potential by a 
set of values on a regular grid = V^(^), from which all 
predicted obser yables are then evaluated by finite differ- 
encing (see e.g. lAbramowitz fc Stegunlll97"2l ). For exam- 
ple, the scaled surface mass density hc is related to ip via 
the Poisson equation, 2k, = V'^t/j (where the physical sur- 
face mass density is S = Scrit and Scrit depends upon 
the angular diameter distances between the observer, the 



lens, and the source). Similarly, the shear 7 = 71 + 172 
and the deflection angle a are derivatives of the potential: 
71 O.5('0,ii - ?/^,22); 72 V^,i2; and a = VV^. 

In our improved method presented here, the regular- 
ity of the grid is relaxed. The advantage of working on 
an irregular and adaptive grid is that it is very flexible, 
which is crucial when dealing with merging clusters. As 
a result we cannot use the standa r d finit e differencing 
formulae of lAbramowitz fc StegunI (|1972l ) , but instead 
extend them to allow for non-u niform differences in the 
spatial variables (see also iDeb et al. 2008 for an alter- 
native approach). We write down the Taylor expansion 

around the point 0^ where we want to estimate the lens- 
ing quantities, in terms of N neighbouring grid point 
locations 0^ (note that 0^ does not have to be one of the 
grid points). This gives 

ni=0 n2=0 1 ^ 

where the maximum derivative order nmax is determined 
by the measurements available to us. For example, if 
our me asurements were to include reduced flexion (e.g. 
iBacon e t al. 2006, Schneider fc Er 2008), then we would 
set rimax = 3 since we would need third order derivatives 
of iIj. 

This gives us a set of N linear equations with M = 
(nmax+iKnmax+2) ^j^j^j^Q^j^g (whlch arc the and its 

derivatives at 0^). We solve this system, requiring N > 
M and using singular value decomposition (SVD) to in- 
vert the matrix. This method is very robust in finding 
a solution for any of the neighbour sub-grid shapes we 
investigated. Note, however, that choosing N = M pro- 
vides rather noisy estimates, especially if nmax is large, 
and occasionally the resulting set of equations cannot be 
inverted. In practice, this situation is easily avoided by 
setting N > M. We repeat the SVD procedure for all 
points 0^ where we either have weak or strong lensing 
galaxies, since this is where we need to be able to predict 
deflection angles and reduced shear ^ = 7/(1 — /^). The 
resulting solution for the spatial derivatives of ip is used 
to evaluate a^, 7, and a. We do not measure higher or- 
der distortions such as reduced flexion from the current 
data, although the method described above can clearly 
be applied to those quantities as well. Note that the SVD 
procedure is only run once, to derive the numerical for- 
mulae for the predicted observables as a function of the 
gridded potential. All subsequent iterations where the 
potential values are varied make use of these formulae. 

The independent strong and weak lensing data are then 
combined by multiplying their likelihoods together. We 
compute the difference between the data (the po- 
sitions of the multiple images, and the ellipticities of 
the weakly lensed galaxies) and their model predictions. 
Since the weak lensing data is noisy and we choose to 
work on a many-parameter lens potential pixel grid, reg- 
ularisation needs to be employed. This process ensures 
that unphysical pixel-to-pixel variations in surface mass 
density are suppressed, and corresponds to asserting a 
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smoothness prior on the reconstruction. In practice we 
compare the current surface mass density and shear map 
with those of an initial model, and penalize strong devia- 
tions in their mean squ are diff erence by adding this term 
to the data (see Bradac et al. 2005 for details). The 
regularisation is chosen (and tested to be) such as to give 
negligible bias in the resulting aperture mass estimates. 

The final global merit function (including data and 
regularisation term) is equal (up to a normalising con- 
stant) to the negative logarithm of the posterior proba- 
bility density function (PDF) for the model parameters 
P, the gridded potential values: 



logP: 



Xwl(V^) , Xsl(V^) 



r]R{ilj) + constant (2) 



(see lBradac et al l|2005l for details). The posterior peak 
values of the potential i/j^ are found by solving the non- 
linear equation dlogP/dt/jk = 0. We linearize this set 
of equations and reach a solution in an iterative fash- 
ion (keeping the non-linear terms fixed at each itera- 
tion step). This requires an initial guess for the grav- 
itational potential; the systematic effects arising from 
various choic e s of this initial model were discussed in 
iBradac et aP (|2006[ ). 

We do not evaluate the covariance matrix of the mul- 
tivariate Gaussian, arguing that the mass reconstruction 
is actually limited by systematic effects rather than sta- 
tistical uncertainty. One could numerically evaluate the 
second derivatives of log P at the peak, or perhaps esti- 
mate the errors on the final reconstruction usi ng MCMC 
sampling ( Marshanl [2QQ6l . lJullo fc Kneibll2009D . However 
with A/grid ^ 3000 this is computationally impractical at 
the moment and is beyond the scope of this paper. 

Instead we focus on systematic uncertainties in the 
magnification in the inner (3^x3^) region around the main 
cluster, where we later search for high-redshift galaxies, 
and which are dominated by unmodeled substructure. 
We investigate the errors on magnification in Sect. [5] by 
adding small scale mass clumps (as well as varying the 
smooth mass profile) to the reconstruction and investi- 
gating their effect. This is adequate for this work: the 
accuracy with which our lens potential model predicts 
the positions of multiple images gives us an indication of 
the range of allowable perturbations. 

The choice of particular grid geometry, the regularisa- 
tion parameter, and the hyper-parameters that set the 
relative weighting between the contributions to 
become critical when weak lensing data on large scales 
(> IMpc) are included and we need a full-field mass re- 
construction. This is not the case in this work, as we are 
only interested in the magnification of the inner region. 
We leave this investigation to future work, where a com- 
plete mass reconstruction of the cluster will be presented. 

3. OBSERVATIONS AND DATA REDUCTION FOR THE 
BULLET CLUSTER DATA 

ACS/WFC imaging of the Bullet Cluster was carried 
out in Cycle 13 (proposal 10200, PI Jones) on 2004 Oc- 
tober 21 and Cycle 15 (proposal 10863, PI Gonzalez) on 
2006 October 12-13. The data consists of two pointings 
centered on the main cluster and the subcluster with 
a small overlap between them. The main cluster has 
been imaged in F606W (hereafter V, 2340s), F775W (i, 
10150s), and F850LP (z, 12700s), while the subcluster 



has been imaged in F435W(2420s), F606W (2340s), and 
F814W (7280s). For the high redshift galaxy survey de- 
scribed below we primarily use the main cluster's Viz 
data. The rest of the data was used in addition when 
generating t he weak and s t rong lensing data catalogs as 

described in IBradac et a l.l (12006^ 

We use the Multidrizz le iKoekemoer et al.l l2002l ) 
routine to align and combine the images. To register 
the images with the astrometric accuracy needed for 
lensing analysis, we determine the offsets among the 
images by extracting high S/N objects in the individ- 
ual, distortion corr ected exposures. We use SExtractor 
(|Bertin fc ArnoutsI 11996 ) and the IRAF routine geomap 
to identify the objects and calculate the residual shifts 
and rotation of individual exposures, which were then fed 
back into Multidrizzle. We use "square" as the final 
drizzling kernel and an output pixel scale of 0.03 arcsec; 
this is smaller than the original pixel scale of the ACS 
CCD, allowing us to exploit the dithering of the observa- 
tions and improve the sampling of the PSF. The limiting 
magnitudes were estimated using 0.2'' x 0.2'' appertures. 
The 5 — <j limiting magnitudes are 28.2 for V, 27.6 for i 
and 27.4 for z-band data. 

4. IMPROVED MASS RECONSTRUCTION FOR THE 
BULLET CLUSTER 

For the strong and weak lensing mass reconstruction, 
we use the deep HS T/ ACS images described above and 
the data used in Bradac et al. (2006) and Clowe e t al] 
( 2006 a ). We use the cat alog of strongly lensed images 
from [Bradac et all (|2006f ), with the addition of two new 
multiple imaged systems and one additional image (see 
Table [J). The we ak lensing catalogs are taken from 
I Clowe et al.l (|2006a) and are obtained from a combina- 
tion of ACS and ground-based data and extends to a 
large FOV; we use the inner 9' x 9' here. The positions 
and redshifts of the strongly-lensed images are given in 
Table [1] (see also Figure [1]). The first new multiply im- 
aged system (Table [H K ) is th e IRAC-selected syste m re- 
ported in Bradac et al.l (120061 ) for which Gonzalez et al.l 
(2009) estimated a photometric redshift of z = 2.7 ± 0.2 
and also found an additional image. The second new 
system (Table [TJ L) is a multiply imaged pair at z ~ 5 
discovered in our search for 2; ~ 5 — 6 galaxies and de- 
scribed in Section[5l We also include a new counter image 
for system F, which we found in close proximity to the 
position predicted by our lens model. 

The reconstruction is set up on an initial grid of 15 x 15 
pixels for a field of 9' x 9'. We then adapt the grid as 
follows. We have the highest S/N ratio in areas close to 
the observed multiple images; in addition, the S/N ratio 
declines as we go towards the outskirts of the cluster. 
We therefore refine the circular regions of 1' diameter 
around the centers of the main and the subcluster by a 
factor of 8 (i.e. each grid is split into 8x8 grid points), 
by a factor of 4 in an annuli between 1' and 1.5', and 
by a factor of 2 in an annuli between 1.5' and 2'. We 
also refine cells surrounding 0.5' around each multiply 
imaged system by 16. The resulting grid structure is 
visible in Fig. [2j the smallest cells have sizes of 2.25", 
corresponding to ~ 10 kpc and to the Einstein diameter 
of massive ellipticals. 

Only for system A has a spectroscopic redshift been 
obtained (Mehlert et al.. .2001) . Hence, we first per- 
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Fig. 1. — Color composites of the three multiple image systems: 
F (top) and the two of K (middle) and L (bottom). Viz color 
composite is used for systems F and L. System K is not detected in 
any of the optical bands; instead, we show Spitzer/IRAC 3.6/im- 
4.5/im-8/im images (see also Gonzalez et al. (2009). North is up 
and East is left in each panel. The cutouts are 10^^ x 10" . 



turbations to the model in Section [5l As the two are 
degenerate, the perturbations that produce the shifts in 
the predicted positions of the multiple images larger than 
the rms residuals quoted above will give us conservative 
errors on the model. 
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Fig. 2. — The reconstructed surface mass density k of the cluster 
for a fiducial source at infinite redshift, ^ oo. 



form the init i al rec onstruction using redshifts given in 
iBradac et aP (|2QQ6l ). which are a combination of photo- 
metric redshifts and lensing predictions from a smooth, 
simply-parametrized model. Then, we project all the 
images for each source to the source plane, take the av- 
erage source position and calculate the rms difference be- 
tween the measured image positions and those predicted 
by the model for the average source. We then vary the 
source redshift (while keeping the mass reconstruction 
fixed) around its predicted value given in Bradac et al. 
^2006). The redshift with the smallest rms for each of the 
sources is given in Table [H The final surface mass den- 
sity maps is given in Figure [3l together with the X-ray 
brightness contours. The average rms in predicted im- 
age positions per image (and source) is lA" . The recon- 
struction has improved significantly with respect to our 
earlier map (in Bradac et al. 2006 the average rms was 
4''). The rms residuals depend upon the final adopted 
grid size. Whereas the final grid size could be further 
decreased, we chose not to do so, as the lack of redshift 
information limits our ability to reconstruct these small 
scales. We can achieve similar fits by adding small scale 
substructure close to the images, or changing their red- 
shifts. Hence, in order to improve the mass reconstruc- 
tion even further, we will need spectroscopic redshifts for 
more than one system, which are not presently available. 

The statistical uncertainties in the redshifts of the mul- 
tiply imaged sources can be estimated while keeping the 
mass reconstruction fixed. However, these are small and 
do not represent the true uncertainties. We would need 
to simultaneously vary the lens model and the redshifts 
to obtain true errors; this, however, is extremely com- 
putationally intensive. At present, however we mimic 
redshift uncertainties by adding small scale mass per- 



5. BULLET CLUSTER AS A COSMIC TELESCOPE: LYMAN 
BREAK GALAXIES AT Z - 5 - 6 

In this section we present a search for V and i-band 
dropout galaxies behind the Bullet Cluster. 

5.1. Dropout selection and photometry 

Galaxies were detected in the HST/ ACS Viz imag- 
ing data described in Section E[ we followed a p roce- 
dure almo s t iden tical to that of iStark et al.l (|2009l ) and 
iCoe et ahl (j2006f ). Faint objects were detected from 

a combined Viz image, and colors measured in aper- 
// 

tures of 0.6 diameter. The fo reground galac t ic ex - 
tinct ion was c orrected following iSchlegel et al.l (|1998l ) 
and Cardelli e t al.l (|198 9). We also applied an aper- 
ture correct ion fo r the F850LP image of —0.05 following 
ICoe et all (j2006[ ). The total magnitudes in each band 
were computed using a combination of aperture color 
and SExtractor parameter MAG_AUTO from the red- 
dest filter F850LP. We then selected filter dropouts with 
the same criteria a s used in Stark et al. (2009) (and also 
in Beckwith et al.l (|2006l) ). For V-dropouts we require 
all of the following: 

{V-i> 1.47 + 0.89(i - z)) or {V - i > 2) (3) 

V-i>1.2 (4) 
i-z <1.3 (5) 
S/N{z) > 5 (6) 

and since we do not h ave ^-band data, w e omit the last 
selection criterion from lStark et al.l (|2009l ). This does not 
significantly affect our results, as there are very few pos- 
sible objects that would drop out in V and be detected 
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Fig. 3. — The Viz (main cluster) and F435W-y-F814W (sub cluster) color composite of the 1E0657-56. The ACS images are inset on 
the colour image from NASA Press release 06-297. Overlaid in red contours is the surface mass density k, from the combined weak and 
strong lensing mass reconstruction. The contour levels are linearly spaced with Ahi = 0.1, starting at k, = 0.7, for a fiducial source at a 
redshift of Zs oo. The X-ray brightness contours from the 500 ks Chandra ACIS-I observations ( Markevit ch et an [2006) are overlaid in 
yellow. North is up and East is left, the field is 7 A' x 5.3^ which corresp onds to 2000 x 1400 kpc^ at the redshift of the cluster. The color 
composite was created following the algorithm from lLupton et al.l (|2004l ) and using the image from NASA Press release 06-297. 



in B (e.g. an AGN in a dusty g alaxy). In addition, 
the study of IBouwens et all (l2QQ7h . which we also com- 
pare to, does not use B-band non-detection either. For 
z-dropouts we use: 

i-z> 1.3 (7) 
S/N{z) > 5 (8) 
(^/A^(^606) < 2) or {V-z> 2.8) . (9) 

In addition, we rejected stars using the IMCAT software. 
We determined the significance and size of each object 
in each passband by convolving the images with a series 
of Mexican hat filters and determining the smoothing 
radius, rg, at which the filtered objects achieved maxi- 
mum significance. This pro cedure, usua l ly u sed for weak 
lensing analyses (see e.g. iClowe et al.l l2QQ6b). gives a 
better estimate of the size than SExtractor. We then 
identified stars from the rg vs. magnitude diagram and 
rejected the bright stars {rriz < 22). Finally, we visually 
examined each object from the catalog in the images, 
removing any additional faint stars, artifacts (at the de- 
tector edges) from the sample. The visual inspection is 
successful in removing stars with {niz < 26), at mag- 
nitudes 26-27, however, t he fraction of stellar co ntam- 
inants is negligible (Bouwens fc Illingwor thl [2006h . Our 
final sample comprises 20 V and 4 i-band dropouts, their 
positions are given in Table [2] and the cutout images are 
shown in Figures [4] and [H 

5.2. Estimate of the completeness and ''lensed" number 

counts 



To estimate the completeness of the sample we ran 
simulations using IRAF task art data. We generated 
1000 objects each at the average z-band magnitude rriz 
of {24, 25, 26, 27, 28}, with colors corresponding to a tem- 
plate spectrum from Bruzual A. & Chariot (1993) for a 
low metallicity {Z = OAZq) starburst galaxy with an age 
of 125 Myr and no dust. We assigned random redshifts 
to the galaxies (uniformly distributed in 4.0 < z < 5.9 
for V, and 5.0 < z < 7.6 fo r z-band dropouts), and then 
follow Madau et al.) (|1996l ) in applying a corresponding 
IGM absorption correction. The galaxies were added to 
the Viz images of the cluster and detected following the 
same procedure. The completeness levels are > 95% for 
ruz = 24, 90% for = 25, 50% for = 26, and 5% 
at rriz = 27. To estimate the completeness for each indi- 
vidual object having rriz < 27, a fourth order polynomial 
as a function of rriz was fitted through above estimated 
completeness levels. The objects were then later added 
to the appropriate bins, assuming these completeness fac- 
tors. We exclude any objects having (lensed, measured) 
magnitudes rriz > 27 when calculating number counts of 
objects as a function of magnitude (see Fig. [6]). The size 
of the field is 10.24 arcmin^. 

We can now calculate the surface densities for our two 
dropout samples and compare them to values from the 
literature (we choose iStark et al.l l2009l . IBouwens et al.l 
120071 . iBeckwith eFall 120061 as thev use the same filter- 
set and a comparison is straightforward). We bin our 
data in two bins of rriz = [25 — 26] and [26 — 27] for 
V and one bin rriz = [25 — 26.5] for z-dropouts. For 
the comparison with the literature the bins are then di- 



TABLE 1 

The properties of the multiply-imaged 
systems used in this work. 





Ra 


Dec 


^pred (^pred06 ) 




104.65852 


-55.950585 




K 


1 04 fi^447 


— 55 Q51 Qn4 






104.63917 


-55.958060 




104.65024 


-55.961577 




L 


104.64327 


-55.963856 


5.7 (N/A) 




104.64694 


-55.958465 






104.65209 


-55.956241 


0.9 (0.8) 




104.66497 


-55.951324 




104.63316 


-55.941395 




A 




C C Q/l Q7QO 


3.24 (3.24) ^ 




104.62954 


-55.941844 




B 


104.63042 


-55.941474 


3.9 (4.8) 




104.63775 


-55.941851 




C 


104.63338 


-55.945324 


2.3 (2.1) 




104.64709 


-55.943575 




D 


104.63528 


-55.951836 


1.5 (1.4) 




104.64008 


-55.950620 




E 


104.64232 


-55.948784 


0.9 (1.0) 




104.56568 


-55.939832 




G 


104.56402 


-55.942113 


1.1 (1.3) 




104.56417 


-55.944131 




104.56293 


-55.939764 




H 


104.56133 


-55.942430 


2.1 (1.9) 




104.56189 


-55.947724 




104.56186 


-55.946114 




I 


104.56052 


-55.942930 


2.3 (2.1) 




104.56141 


-55.944264 




104.56909 


-55.946016 




J 


104.57025 


-55.944050 


1.3 (1.7) 



^ First two images are from Brad ac et al.l f2006), 
the last was predicted in its close proximity by 
the lens model. 

^ Only image A has a measured spectroscopic 
redshift (Mehlert et al. 2001). Others are the 
photometric redshifts refined by the gravitational 
lensing model. Systems from here on were used 
already in Bradac et al. (2006), the redshifts used 
there are listed in parenthesis. 

vided by 2 and 3 respectively to obtain counts in 0.5- 
mag bins. Our results are presented as grey open circles 
with dashed errorbars in Fig. [6] (the determination of the 
"unlensed", true number counts — shown as black solid 
circles in Fig. [6] — will be described in Section The 
errors include both Poisson errors and sample variance 
(also commonly referred to as cosmic variance). To cal- 
culate the sample varia nce we use the prescription from 
iSomerville et al.' (2004) and the correlation function from 
^agamine et al. (2007). The sample variance for a single 
ACS field at z ~ 6 is 50%, whic h is also in agreement wit h 
sample variance calculator by iTrenti fc StiavellH (|2008l ) . 
We immediately note that the lensing effect has increased 
the number of detected dropouts compared to blank sur- 
veys with similar or better exposure times to ours. Even 
though gravitational lensing effectively makes the actual 
observed solid angle smaller (due to magnification), the 
effective slope of the luminosity function ^ is steep at 
these magnitudes (— d(log <l>)/d(log L) > 1), hence we in- 
crease the expected number density of sources by going 
deeper. 

5.3. Number counts including magnification correction 

To evaluate the true (unlensed) number counts, we 
use the mass reconstruction described in Section [4] to 
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Fig. 4. — Cutouts of V-band dropouts. The cutouts are shown for 
filters Viz (in two columns), labels correspond to labels in Table [2] 
The sizes are 10'' x 10''. 

estimate the effective change in survey area. We use 
the defiection angle map for a redshift z = 5 source 
for the V-dropouts, and for a z = 6 source for the i- 
dropouts. We project the observed field from the lens 
plane to the source plane pixel by pixel and measure the 
corresponding change in solid angle by simple numerical 
integration. This resulted in a fractional area decrease 
of 0.25 ± 0.02 for z = 5 and ^ = 6 case. The total 
shrinking of the field changes little with redshift, as the 
lensing strength (or angular diameter distance ratio be- 
tween lens and source, and observer and source) at these 
redshifts f or a lens at z ^ 0.3 i s nearly constant (see 
e.g Bartelmann fc Schneiderl l2001f ). The uncertainty was 
estimated by evaluating the area change using the mag- 
nification error map (described below). Neither a source 
redshift change (supported by the calculated solid angle 
change for sources at z = 5 and 6) nor the magnifica- 
tion errors have a significant effect on the survey area 
shrinkage. The position of the critical curves changes by 
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Fig. 5. — Same as Fig. |4l except for i-band dropouts. 

~ 5 arcsec between z = 5 and z ^ oo, which is larger 
than the expected accuracy; even so, this has an almost 
negligible effect (within the errors quoted above) on the 
total solid angle change of the field. 

Next we need to evaluate what the apparent magni- 
tudes of the sources would have been in the absence of 
lensing. For this, we evaluate individual magnifications 
at the positions of these sources using our lens model. 
The errors on each magnification come from the errors 
on the mass reconstruction, and typically increase as the 
image approaches the critical curve. These errors are 
dominated by small-scale substructure not accounted for 
in the lens modeling, and by erroneous redshifts of mul- 
tiple images used for the reconstruction. They are less 
affected by the typical uncertainties in cluster mass pro- 
file as discussed below. 

5.4. Errors on magnification 

To estimate the errors we used (for convenience) the 
approximate parametrized model for the Bullet Clus- 
ter (consisting of one PIE MD (pseudo-isotherma l ellip- 
tical mass density, see e.g. iLimousin eTaDHQQll) com- 
ponent for each of the main and subcluster, as well as 
30 SIS (Singular Isothermal Spheres) placed over the full 
mass-reconstruction field representing the cluster galax- 
ies (these are the confirmed cluster members in the inner 
5' X 3' field). These galaxies were taken to have line-of- 
sight velocity dispersion (Jmemb following (Jmemb oc L^/^ 
and a fiducial galaxy with a F606W-band magnitude 
^F606W = 18 has crmemb,i8 = 250kms~"^. We chose 
singular rather than non-singular isothermal spheres, as 
they will create larger local magnification changes. In ad- 
dition, strong lensing studies of galaxy-scale lenses show 
that they have steep mass pr ofiles in the centre, and 
hence small core radii (see e.g. lKeeton|[2QQ3f ). 

The above model gives us a manageable smooth model 
for the cluster mass distribution. To model the errors 
on the magnification due to substructure, we also ran- 
domly distributed 50 mass clumps with position covering 
our Viz images (where dropouts were searched for) and 
magnitudes between mpeoew = 22 - 24 (and the same 
scaling of mass as used above). Again, we are trying 
to maximize the magnification perturbations, hence us- 
ing singular profiles is conservative. The mass fraction 



in these substructures was ~ 8%; this is purposely high, 
since we want to be conservative in estimating the magni- 
fication errors. This model does not reproduce the mul- 
tiple image positions in detail (i.e. we did not search for 
the best-fit model - the substructure fraction would have 
been lower and systematic errors due to uncertainties in 
redshifts and the main lens profile would be underesti- 
mated); however this is not needed for the error esti- 
mation. This additional substructure changes the criti- 
cal curve positions by 5 — 10 arcsec, which is more than 
the uncertainty with which multiple images are recov- 
ered (in this and in other well-modeled clusters, see e.g. 
lJullo fc KneibI [2QQ9h . We generated 100 reahsations of 
cluster substructure in this way, and then evaluated the 
magnification error as the standard deviation of the mean 
magnification at each pixel. The map of standard devi- 
ation in the change of magnitude (Jmag — ^^1 (m) as a 
function of position is given in Fig. [71 On average the 
errors are small across the observable field, only in the 
vicinity of the critical curve the errors become large. The 
errors at the positions of the dropouts are given in Ta- 
ble [2 

In addition, the average error ((Jmag) as a function of 
average magnification (/i) is shown in Fig. [H As a com- 
parison we also give the cumulative area as a function 
of magnification /i, this time calculated for our recon- 
structed magnification map from Section [3l These plots 
illustrate that less than 0.1 arcmin^ of the survey (i.e. 
< 1% of the area for a single pointing ACS observation) 
has magnifications /i > 10 and on average low errors. 
It is true, that we are more likely to find images in re- 
gions of large magnification; however only 3 out of 20 
systems have highly uncertain magnifications (see Ta- 
ble [2]). Hence the errors on the positions of the magni- 
tude bins in evaluating the number counts (Fig. [6|) are 
much smaller than the bin widths. 

Above investigations use a fixed smooth model to de- 
termine uncertainties. To further evaluate uncertainties 
due to profile uncertainties, we have investigated magni- 
fication properties of three very different smooth profiles, 
singular isothermal sphere (SIS), non-singular isothermal 
sphere (NIS) and Navarro Frenk and White (NFW) pro- 
file. NFW and NIS profiles give rise to two critical curves, 
radial and tangential, whereas SIS's radial critical curve 
shrinks to zero size. The existence of radial arcs and 
simulations tell us that SIS is not a valid description for 
cluster mass distribution. Systems with > 2 multiply 
imaged systems with known redshifts (like the one used 
here), and even more so, with included weak lensing in- 
formation c an distingui sh between thes e profiles (see e.g. 
Sand etalj 120081 . [Limousin et al.' 2008) . 

Still, we have investigated their magnification prop- 
erties. First we match the Einstein radii of all three 
models, which is easily measured in practice even with 
a single multiply imaged system with known redshift. 
If also matching the radial critical curve, the change in 
area (due to magnification) is < 1%. If only the tangen- 
tial critical curve is matched, and the radial ones differ 
by ~ 10 arcsec these changes are still below < 20%, only 
when we compare SIS with NFW is the change in area 
significant. It is only significant for small observing fields 
(like the one used here) and can be as high as 50%, how- 
ever this case is unrealistic as noted before. Furthermore, 
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for a field twice tiiis size, tiie errors due to uncertainties 
in smooth model are again below 5%. The individual 
magnifications differ by < 20% everywhere outside the 
tangential critical curve and are mostly in 5% regime, 
which is where most of our sources are located. 

In conclusion, the errors due to lensing are smaller than 
the sample variance. Taking into account sub-dominant 
errors, the increased number of objects detected, and the 
ability to reach otherwise unobservably faint magnitudes, 
it is advantageous to use well-modeled clusters as cos- 
mic telescopes to measure the average properties of back- 
ground {z > 5) objects, including their luminosity func- 
tion. This conclusion is not limited only to high-redshift 
sources. The system K is at z ~ 2.7, but magnification 
factors of /i ^ 50 make this ordinary (i.e. not ultra- 
luminous) sub-mm so urce detectable with Spitzer/IRAC 
([Gonzalez et al."2009'), AzTEC (At acama Submillimeter 
Telescope Experiment, Wilson et al. 2008 and BLAST 
(Balloon-borne Large- Aperture Submillimeter Telescope, 
iRex et al.ll2009l ). 

6. CONCLUSIONS AND OUTLOOK 

We have presented a new mass reconstruction for the 
post-collision Bullet Cluster. In addition, we have used 
its extraordinary capabilities as a cosmic telescope to 
study z ^ 5 — 6 objects. Our main conclusions are as 
follows: 

1. Our new strong and weak lensing mass reconstruc- 
tion method, that uses a non-uniform adaptive 
pixel grid, performs very well in reconstructing the 
cluster potential in the vicinity of multiple images. 
In particular, we have improved the accuracy with 
which the image positions are reproduced to an av- 
erage rms residual of 1.4^'. The limitation in push- 
ing this number to the observed uncertainties in im- 
age position (~ 0.1'') is mostly in the lack of secure 
redshifts (only one spectroscopic redshift is known 
for the 12 multiple image systems used here). The 
method can resolve the small scale substructure 
and reconstruct the images even more precisely, 
once the redshifts become available. 

2. Bullet Cluster is an efficient cosmic telescope. We 
have used its power to search for z ~ 5 — 6 
objects using the dropout technique. After cor- 
recting for the lensing effect and taking into ac- 
count the errors from the imperfect magnifica- 
tion map, we derive the surface densities of V 
and z-band dropouts. They match well the den- 
sities (in absolute value) d erived from previous 
studies (IStark et all 120091 . iBouwens etaU 120071 . 
iBeckwith et al.l l2_^ using either deeper or larger 
surveys (HUDF and GOODS). E.g. GOODS (re- 
lease 1) data has similar exposure times (5000s in 
V and i and 10660 in z) over 360 ar cmin^ field. 
Scaling total numbers of dropouts from lStark et al.l 
(|2009) (using area ratio between our and their 
survey) we would obtain 14 V and 3.5 i-band 
dropouts; whereas we observe 20 and 4. In ad- 
dition they are at fainter (intrinsic) magnitudes as 
their sample. 



3. Our sample numbers are much smaller (as our 
search area is 10 arcmin^ , comp ar ed to the area 
used in e.g. IStark et al.l (12009*), "B ouwens et al.l 
(|2007f ) with > 300 arcmin^), and so the error- 
bars are dominated by sample variance and Pois- 
son noise. However, our results show, that at the 
same observed magnitudes, we gain in the number 
of sources observed compared to a blank field of 
the same size. In addition, we observe intrinsically 
fainter sources than would otherwise be possible in 
a blank field observed to similar depths. Finally, 
multiply imaged sources are easily discriminated 
from contaminants, as lens model roughly predicts 
their redshifts, and hence can be distinguished from 
z ^ 2 dusty galaxy and cool stars. 

This analysis can and will be extended to higher red- 
shift objects (z and J-dropouts). The main requirement 
for such studies is deep optical and near-IR data, as well 
as high resolution magnification maps such as the one 
presented in this paper. These studies can also be readily 
expanded to a larger sample of clusters. With a sample 
of ~ 20 massive clusters, with well understood magnifica- 
tion properties and observed to similar depths as the Bul- 
let Cluster we can increase the expected number counts 
and reduce the errors due to sample va riance and Poisso n 
sampling to those of GOODS (see e.g. IStark et al.ll2009[ ). 
However, due to lensing, we can push observations to at 
least a magnitude deeper than GOODS and only slightly 
shallower than HUDF (in much shorter observing time). 
Finally, highly magnified images (/i > 50) potentially al- 
low for (otherwise prohibitive) spectroscopic follow up 
of these high redshift galaxies. Highly magnified objects 
will often sample populations that would otherwise be 
unobservable by any practical means. 
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Fig. 6. — Surface densities (number per arcmin^ per 0.5 mag) of V (left) and i-band (right) dropouts (circles). Due to small number 
statistics, we bin the data in two bins of ruz = [25 — 26] and [26 — 27] for V and one bin ruz = [25 — 26.5] for i-dropouts and recalculate them 
to 0.5-mag bins, which is the binning used in the literature (crosses). Only sources for which our completeness is greater than 5% (i.e. those 
with ruz < 27) are presented and completeness corrections are discussed in the text. Direct measurements (gray open circles with dashed 
errorbars) have been corrected for the effect of gravitational lensing, giving the intrinsic (unlensed) number counts (black full circles with 
solid errorbars). The errors for our data include sample variance, Pois son errors , and errors due to lensing (the latter of which are small in 
comparison). The crosses are all calculated for the bins indicated by St ark et"aD (2009) data, the other two are shifted by 0.05 and 0.1 mag 
for clarity. Bouwens et al. (2007), and Stark et al. (2009) data have been corrected for comple teness, whereas Beckwit h et al.i (j2006) data 
were not. The Stark et al. (2009) points were calculated using GOODS data, Bec kwith et al.l (2006) only include HUDF number counts, 
whereas .Bouwens et al.. (.2007) use data from HUDF, HUDF-Ps, and GOODS. 




Fig. 7. — The error on intrinsic (unlensed) source magnitude 
(Jmag (grey-shade) as a function of position. Most of the lens-plane 
area has errors < 0.1 mag. The field we study here is approximately 
centered on the main cluster and is given as a box. The map 
assumes a source at z ^ oo, the map for z = 5 or 6 would differ 
very little from the one presented here. 



Barrena, R., Biviano, A., Ramella, M., Falco, E. E., & Seitz, S. 
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Fig. 8. — Top: the average error on magnitude {(Jmag) as a func- 
tion of magnification. For magnifications oi fj, < 10, the errors are 
< 0.1 mag on average. The exact error as a function of location is 
given in Fig. [Jl Bottom: cumulative number of image plane pixels 
having magnifications larger than a certain value. < 1% of the 
pixels will have > 10, which is the only regime where the errors 
are of any significance for determining the intrinsic luminosity of a 
sample. Both plots assume a source at z ^ oo. 
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TABLE 2 

The properties of F-dropouts (top) and z-dropouts (bottom). 





Ra 


Dec 


mz ± Aphot 


(V-i) 


{i-z) 








1 


104.59377 


-55.970039 


26.74 =b 0.07 


2.2 ± OA 


0.69 ± 0. 


,14 


2.68 ± 0.08 


27.8llg:°| 


2 


104.60201 


-55.954700 


27.89 ±0.10 


1.8 ±0.2 


-2.42 ±0. 


,15 


3.97 ±0.26 


9Q QQ+0.12 


3 


104.60582 


-55.969944 


26.45 ± 0.06 


> 1.5 


0.20 ± 0. 


,11 


3.17 ±0.09 


27.701°:°? 


4 


104.61459 


-55.951771 


27.42 dz 0.05 


2.2 ±0.1 


-0.54 ± 0. 


,10 


4.0 ±0.9 


28.93+lil 


5 


104.61855 


-55.940559 


26.24 dz 0.05 


3.3 ±0.5 


0.86 ± 0. 


,13 


5.8 ± 1.2 


28.15_Q 25 


6 


104.63613 


-55.925079 


25.53 lb 0.05 


2.3 ±0.6 


0.59 ±0. 


,10 


3.78 ±0.15 


26.98l°i°? 


7 


104.64101 


-55.935631 


26.97 ±0.05 


2.3 ±0.7 


1.17 ±0. 


,13 




291^ 


8 


104.64129 


-55.964703 


26.38 ± 0.06 


2.3 ±0.2 


0.14 ±0. 


,11 


12 ±2 


29.09l°i? 


9 


104.64133 


-55.966335 


25.04 ± 0.05 


2.5 ±0.4 


0.97 ±0. 


,11 


10 ± 1 




10 


104.64437 


-55.973896 


25.81 ± 0.05 


2.4 ±0.6 


0.87 ±0. 


,10 


5.2 ±0.3 


27.64l°:°l 


11 


104.64603 


-55.969257 


25.68 ± 0.05 


2.8 ±0.3 


-0.11 ±0. 


,10 


6.5 ±0.8 


27.70t°o:\t 


12 


104.65024 


-55.961578 


25.52 ±0.05 


5.7 ±4.0 


0.88 ±0. 


,10 


480 ± 30 


32.22+°°! 


13 


104.65026 


-55.966267 


25.53 ± 0.05 


2.6 ±0.2 


0.18 ±0. 


,10 




2911 


14 


104.65306 


-55.930065 


26.30 ± 0.05 


2.1 ±0.1 


0.12 ±0. 


,10 


4.6 ±0.2 


27.961°;°? 


15 


104.65884 


-55.948055 


27.23 ± 0.06 


4.0 ±0.2 


0.06 ± 0. 


,11 


207 ±3 


OQ riQ + 0.06 

OO. UZ_QQg 


16 


104.66796 


-55.949810 


25.60 ± 0.05 


3.2 ±0.6 


0.72 ± 0. 


,10 




291^ 


17 


104.67584 


-55.945240 


26.01 ±0.05 


3 ±1.0 


0.32 ±0. 


,10 


7± 1.0 


28.19+0 17 


18 


104.67677 


-55.943958 


28.18 ±0.05 


1.5 ±0.10 


-1.42 ±0. 


,10 


7.1 ±0.5 


30.311°:°^ 


19 


104.68509 


-55.951725 


27.14 ±0.05 


3.9 ±0.2 


-0.23 ±0. 


,10 


8.6 ±0.4 


29.471°;°? 


20 


104.68786 


-55.943943 


27.36 ± 0.09 


2.2 ±0.5 


-0.05 ± 0. 


,14 


4.8 ±0.2 


29.06l°;i° 


1 


104.63920 


-55.925922 


26.29 ± 0.05 


> 0.5 


1.33 ±0.10 


4.9 ±0.2 


28.011°;°? 


2 


104.65271 


-55.930984 


25.44 ± 0.05 


> 1.4 


1.33 ±0.10 


4.6 ±0.2 




3 


104.66599 


-55.931286 


25.38 ± 0.05 


2±1 


1.34 ±0.11 


3.4 ±0.2 


26 71+0-08 


4 


104.67469 


-55.929668 


26.04 ± 0.05 


> 0.5 


1.71 ±0.13 


3.1 ±0.1 


z/ .zu_Q 06 



Note. — For the final error on rriz^int (Ai+p) photometric and magnification errors were added in 
quadrature. 
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